Second order perturbations of rotating bodies in 
equilibrium; the exterior vacuum problem. 
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We study the exterior vacuum problem for first and second order stationary and axially 
symmetric perturbations of static bodies. The boundary conditions and their compatibility 
for the existence of an asymptotically flat exterior solution are discussed. 

> 

. 1 Introduction 

O 

Finding global models for rotating objects in general relativity has proven to be extremely 
ly-j ' difficult, even for axially symmetric configurations in equilibrium. So far there are no known 

explicit global models except for spherical stars (hence non-rotating) or the Neugebauer and 
Meinel disc of dust [1], where the matter source is encoded as suitable jumps on the first 
q-i derivatives of the metric, which is otherwise vacuum everywhere. 

^ 1 Since a proper understanding of rotating objects in equilibrium within the context of 

(^i), general relativity is obviously fundamental for many astrophysical situations, alternative 

methods have been developed over the years. The two most important ones are, without 
doubt, the use of numerical methods and perturbation theory. The former can handle fully 
relativistic situations with intense gravitational fields and high velocities and has produced 
a large variety of very interesting results. However, the fact that closed expressions are never 
found in this setting leaves plenty of room for other approaches, like for instance perturbation 
theory. Its field of applicability is restricted to slowly rotating stars (so that the perturbation 
parameter can be taken to be the maximum angular velocity of the star, for instance) and 
also a wealth of results have been obtained (see e.g. [2] and references therein). 

This short note is a summary of a longer paper [3] by the same authors where many more 
results and their proofs can be found. The object of this contribution is to explain in a few 
words the main motivation of this work, present the main result and give some indications 
of how it can be proven by analyzing the simplest possible situation. 

The aim of this work is to study perturbation theory of rotating stars from a different 
perspective than normally done. We want to consider slowly rotating bodies with an arbi- 
trary matter content and we wish to concentrate on the effect that the rotation has in the 
vacuum outside region. Thus, our background spacetime is composed of a static object with 
some non- vanishing energy- momentum tensor (typically a perfect fluid, but not necessarily). 
Although we have in mind the case when the static body is spherically symmetric, which 
is the physically most relevant one, all our results hold also for axially symmetric back- 
grounds. We assume, as it is usual in this context, that the object has a sharp boundary and 
that there is no matter outside it, so that the metric exterior to the body is vacuum. The 
metrics inside and outside the body must satisfy certain junction conditions (see [4] for a 
detailed account) in order to produce a well-defined spacetime. Given this background, we 
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want to perturb it arbitrarily in the interior except for the restriction that the object is still 
in equilibrium and has axial symmetry (i.e. the perturbed metric is stationary and axially 
symmetric). Furthermore we want to do perturbations up to second order in perturbation 
theory. The necessity of going to second order comes from the fundamental results obtained 
in the seminal paper by Hartle [5] where rigidly rotating perturbations of spherically sym- 
metric and static perfect fluids were analyzed and where it was found that to first order 
rotation only affects the , (j> crossed term in the metric (more technically, it only produces 
so-called axial or odd perturbations) without modifying, for instance, the spherical shape of 
the object, while second order perturbations, already affect the shape of the body and the 
rest of the metric components. 

In Hartle's paper some heuristic arguments were used at some places, specially regarding 
the matching procedure of the perturbed metrics. One of the aims of our paper is to set 
up a proper theoretical framework so that, in a future work, we can make rigorous all the 
arguments used by Hartle. This is important insofar Hartle's paper has served as the basis 
of many developments in perturbation theory of rotating objects over the years. 

2 Brief summary of stationary and axially symmetric rotating 
bodies 

A spacetime describing a stationary and axially symmetric rotating body with a boundary 
surface is composed by two regions; one inside the body solving the Einstein field equations 
with matter, and another outside the body solving the vacuum field equations and being 
asymptotically flat (because we are dealing with an isolated body). Furthermore the two 
metrics must fulfill the so-called matching conditions on the boundary of the body, i.e. on 
a timelike, stationary and axially symmetric hypersurface S. The vacuum field equations 
outside the body can be written in terms of a coupled system of non-linear elliptic PDE for 
two scalars U and Q called the "norm" and the "twist" potentials respectively and which are 
defined in terms of the unique stationary Killing vector £ which is unit at infinity. The twist 
potential is defined only up to an arbitrary additive constant which is fixed by demanding 
fl — > at infinity. With this choice, J7 vanishes if and only if the spacetime is static, hence 
the twist potential determines whether the body is rotating or not. In the vacuum region, 
there exist local coordinates {t, <f>, p, z} called Weyl coordinates which are adapted to the 
stationary and axial Killing vectors and which locate the axis of symmetry at p = 0. This 
coordinate system is defined uniquely except for an additive constant in z, which can be 
fixed whenever the spacetime has an equatorial plane by choosing z — on the equator. We 
will assume that the Weyl coordinates exist also globally in the vacuum exterior region (see 
[6] for global existence results of the Weyl coordinate system) . In this setting U and fi are 
functions of p and z alone and satisfy the so-called Ernst equations, which read 

A 7 [/+i e - 4l/ {d(2,df2), ( = 0, 
A 7 i? - 4 (<M, dC/) 7 =0, 

where A 7 is the Laplacian of the flat metric 7 = dp 2 + dz 2 + p 2 d(j) 2 and ( , ) 7 denotes scalar 
product with respect to it. The asymptotic flatness condition demands U — l—M/r+0(r~ 2 ), 
f2 = —2zJ/r 3 + 0(r~ 3 ), for some constants M, J and where r 2 = p 2 + z 2 . The boundary 
of the rotating body as seen from the exterior region will be denoted by S E and will be 
defined by two functions p = p{p), z = z(p). If the metric inside the body is assumed to be 
known, then the matching conditions can be seen [7] (see [8] for the complete generalisation) 
to be equivalent to the following data: (i) the explicit form of the matching hypersurface 
in Weyl coordinates, i.e. the functions p(p),z(p) and (ii) the values of U and ft (the latter 
except for an additive constant) together with their normal derivatives on E E . Notice that 
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the Ernst equations are elliptic, which means that appropriate boundary data consist of 
fixing the value of the functions or their normal derivatives on the boundary (or perhaps a 
combination of both), but never the full values of the functions and their normal derivatives. 
In more technical terms the boundary conditions are of Cauchy type, which is unsuitable for 
elliptic problems. This property reflects the fact that given an arbitrary metric describing a 
rotating body in equilibrium, in general there will not exist a vacuum solution matching with 
the given metric and extending to infinity in an asymptotically flat manner. The problem 
of finding a global model of a rotating object is truly global in nature and cannot be broken 
into an interior and exterior problem without paying some price. In our case this translates 
into the necessity of dealing with an overdetermined boundary value problem for an elliptic 
system. Thus, one has to worry about existence of the solutions, i.e. which are the restrictions 
that the boundary data (and hence the interior metric) must satisfy so that they truly 
represent an isolated rotating body. Let us stress here that uniqueness is a much simpler 
problem which can be solved in full generality [7] . With regard to existence, there are results 
involving an infinite set of compatibility conditions on the boundary data which are necessary 
for existence to hold [9] . Whether they are also sufficient is still an open problem. 

3 First and second order perturbations of the exterior region 

After this brief reminder of the non-linear case, let us move into perturbation theory. Calling 
e the perturbation parameter, we consider a one-parameter family of spacetimes depending 
on e. Since we take every element in the family to be a stationary and axially symmetric 
spacetime, vacuum outside some spatially compact region (defining the rotating body) and 
asymptotically flat, we have two families of functions U £ (p, z), L7 e (p, z), which for all e satisfy 
the Ernst equations (1). Defining U' = d e U € \ e= o, U" = <9 e <9 e £/ e | e= o an d similarly for Q' and 
Q" , and using now that the background spacetime is static, i.e. that /2 e | e=0 = 0, we find the 
first and second order perturbed field equations 

A 7 U = \ A 7 I# + e- AUa (df2' , d^) 7 = 

AX - 4 {dfr , d(7 ) 7 = / ' A 7 r# - 8 (df2' , dU^ - 4 (d/2£, dC/ ) 7 = 

simply by taking first and second partial derivatives of (I) with respect to e and evaluating 
the result at e = 0. 

Notice that we are performing the perturbation (i.e. derivative with respect to e) by 
considering p, z as independent of e. This entails a suitable identification between the dif- 
ferent spacetimes for different e. Such an identification must always be made in order to 
define metric perturbations. However, the identification is not unique, as we could perform 
an arbitrary diffeomorphism on each element of the family V e of spacetimes and the diffco- 
morphisms can obviously depend on e. This freedom in the identification implies the gauge 
freedom which is inherent to perturbation theory in general relativity (and indeed in any 
geometric theory). Thus, when we take perturbations by performing derivatives with respect 
to e with fixed {p, z} we are effectively fixing the gauge. 

We should consider now which is the domain where the equations (2) hold and what 
kind of boundary data need to be fulfilled. Given the interior family of metrics, the matching 
conditions fix (for every e) two functions p e (p) and z e {p) defining the matching hypersurfacc 
(i.e. the boundary of the body) for every e. They can also be seen as two-surfaces defined 
in Euclidean 3-space E 3 with the flat metric 7 written in cylindrical coordinates. If e is 
close enough to e = we have a family of axially symmetric surfaces E e in E 3 , all of 
them diffcomorphic to the surface corresponding to the static background, which will be 
denoted by Uq. We will furthermore assume that the background surface is diffcomorphic to 
a sphere, which is physical reasonable (and certainly true whenever the static background 
is spherically symmetric). 
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Let us also choose the range of variation of p, so that there exist two fixed values ps and 
fj> 7v such that p e {ps) = Ph{pn) = for all e, i.e. that the intersection points of the surfaces 
with the axis of symmetry occurs at the same value of p. As explained before, the matching 
conditions together with the interior metrics provide us with four functions, which we denote 
as G e , L e , Y e and W e , all of them defined on S e and such that the boundary values for the 
exterior problems are U € \z € — G € , n e (U e )\x: e = L e , f2\s e — Ye, n e(^e)\s e = W € , where 
n e = —z € d p + p e d z \s e is a vector field orthogonal to S e (dot denotes derivative with respect 
to p). Notice that in these expressions the right-hand sides are functions of p and e alone, 
while for the left hand sides they are functions of E 3 evaluated on a (moving) surface S € . We 
can now take derivatives of these expressions with respect to e (at p constant) in order to 
obtain the boundary values for the perturbed functions Uq, fl' , Uq and £2q and their normal 
derivatives on the unperturbed surface S . In particular, it follows that the domain where 
the perturbed Ernst equations hold coincides exactly with the domain corresponding to the 
static background. Moreover, the boundary data of U', U", Q' and Q" on the unperturbed 
hypersurface Eq can be explicitly calculated in terms of the interior background metrics and 
the interior first and second perturbation metrics alone. The resulting expressions are long 
and will not be given here (see [3] for a detailed discussion and explicit expressions). Having 
discussed briefly the perturbed exterior problems and their boundary conditions, we can 
now address the issue of existence of the exterior solution fulfilling these boundary data. 

4 Compatibility conditions 

In the previous section we have seen that, as in the non-linear case, the perturbed exterior 
problem involves an elliptic system of equations with Cauchy boundary data. Again this is 
an overdetermined system and we should not expect asymptotically flat solutions to exist for 
any Cauchy data (i.e. for any interior perturbation). Thus we need to address the question of 
what is the set of necessary and sufficient conditions that the boundary data must satisfy so 
that solutions exist. Asymptotic flatness demands lmv^oo Uq ~ linv^oo f2' = linir^oo Uq = 
linv^oo J?o = 0- The perturbed Ernst equations can be collectively written as 



where u = u(p,z) stands for Uq, Uq, etc., and j = j{p,z) represents the inhomogeneous 
terms in the second order perturbation equations. The metric 7 corresponds to either 7, 
for the [/-equations, or 7 = e~ 8Ua j, for the i?-equations. The domain (Dq,^) is clearly un- 
bounded because 7 is an asymptotically flat metric. Thus, the compatibility conditions 
for the boundary values of Uq, Uq, Q'q, Qq can be studied as particular cases for the 
compatibility conditions of the Cauchy problem for the general inhomogeneous Poisson 
equation (3) defined on an unbounded asymptotically flat region (.Do, 7) with boundary 

= {p = po(p),z = zo(fx),(j) = (p}, corresponding to the boundary of the static back- 
ground metric. Furthermore, we will assume that j tends to zero at infinity at least like 1/r 4 
(which follows in our case from asymptotic flatness). 

In order to give a flavor of why Theorem 1 holds, let us concentrate on the simplest 
possible case, i.e. when 7 = 7 and the source term j vanishes. A simple consequence of 
Gauss' theorem is the so-called Green identity, which reads: for any compact domain K e Dq 
and any function ip (both suitably diffcrcntiable) 



where n 7 is a unit (with respect to 7) normal vector pointing outside K, rj 7 is the volume 
form of (.Do, 7) and d5 7 is the induced surface element of dK. We intend to apply this 
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identity to a function tp that (i) solves the Laplace equation A 7 ?/> = on Do, (ii) admits a 
C 1 extension to dD n = S and (iii) it decays at infinity in such a way that rtp is a bounded 
function on D . A function tp satisfying these three properties is called a regular 7 -harmonic 
function on Dq (if a function satisfies just (ii) and (iii) and is C 2 on Do we will call it regular). 
For such a function we may take K = Do in (4) because the integral at the boundary "at 
infinity" can be easily shown to vanish. Thus, denoting the boundary data of u on So as 
u\e = fo and n 7 (w)|i; = fi of u, the expression above becomes 



Is, 



V/i-/on 7 (V) dS 7 = 0. (5) 



These are obviously necessary conditions that the overdetermined boundary data must sat- 
isfy in order for a regular solution u of the Laplace equation to exist. It is natural to ask 
whether such conditions are also sufficient. More precisely, give two arbitrary continuous 
functions fo and f\ on So which satisfy (5) for any choice of regular 7-harmonic function tp. 
We want to check whether there always exists a function u satisfying the Laplace equation 
A 7 u = 0, with ru bounded at infinity and such that the boundary equations u\s = fo, 
n j( u )\s — h are satisfied. The answer is yes as we show next. Consider the Dirichlet 
problem A 7 it = with u\s a = fo- Standard elliptic theory tells us that this problem al- 
ways admits a unique solution u which tends to zero at infinity. Let us define fi on S as 
fi = n 7 (u)|£ . Since u solves the Laplace equation, expression (5) must hold replacing fi 
with f\. Furthermore, our assumption is that fo and f\ satisfy (5) for any regular 7-harmonic 
function tp. Subtracting both expressions we get f s V(/i — /i)d<S 7 = 0. However, since the 
regular 7-harmonic function tp is arbitrary and the Dirichlet problem for the Laplace equa- 
tion always admits a solution, we have that ip\s is an arbitrary continuous function. This 
readily implies /1 = f\ and hence compatibility of the overdetermined boundary data. 

However, condition (5) has a serious practical disadvantage, namely that it must be 
checked for an arbitrary decaying solution tp of the Laplace equation. This makes it not 
useful in practical terms. Our aim is to reduce the number of solutions tp that must be 
checked in (5) while still implying compatibility of f and f\. Here is where axial symmetry 
plays an essential role. This allows us to reduce the number of conditions to be checked to 
just a one-parameter family set. Without going into the details, let us just state that, in the 
axially symmetric case, conditions (5) restricted to the set of functions 

V>,*) = l/vV + (2-2/) 2 , (6) 

where y is a constant such that the point {p = 0, z = y,<p} lies outside the domain Do, 
are already necessary and sufficient for existence of the solution u. Thus the number of 
conditions reduces dramatically and becomes manageable. The main result we present in 
this contribution is the generalization of this result to the four cases corresponding to the 
perturbed Ernst equations. The integrals to be performed are, of course, more complicated 
in general but the idea is still the same. In order to write down our main theorem we 
need to introduce some notation. First of all, define z$ < zn as the values of z at the 
intersection of So with the axis of symmetry (i.e. the south and north poles respectively) 
and restrict y to the interval (zs,zn)- We define an angle T y e [0, 2tt) by cosT y (p, z) = 
(z — y)/y/f? + (z — y) 2 , s'mYy(p, z) = p/ \/ p 2 + (z — y) 2 and three functions W y , Q+ and 
Q_ as the unique vanishing at infinity solutions of the following compatible PDE 

dW y — cos Y y dllo + sinT y * dUo, 

dQ± = e - 2Uo±2W y [(=pl - cosr v )dO' - smT^d^] . 



Here * means Hodge dual with respect to the metric dz 2 + dp 2 . Then, we have 
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Theorem 1. Let fo, fi be continuous axially symmetric functions on a C 1 simply connected, 
axially symmetric surface Eq o/E 3 . Let this surface be defined in cylindrical coordinates by 
{p = Po(p), z = z o(p),(f) = <p}, where p, e [^s^Pn] and ps < p^ are the only solutions of 
Po{p) = 0. Call zs = z(ps) and zjq = z(ppj) and assume z$ < zjq (i.e. that these values 
correspond to the "south" and "north" poles of the surface, respectively). Denote by D the 
exterior region of this surface and n = — z d p + pod z a normal vector to it. Then, 

(i) the Cauchy boundary value problem A 7 £/q = 0, Uq\e = fo, n (Uo) \s = A admits a 
regular solution on D if and only if 

/ bPy fi ~ fo n(V'y)] Podp, = 0, Vy e (z s , z N ), 

J ps 

(ii) the Cauchy boundary value problem A 7 I2 — 4 (di? , d?7o) 7 ~ 0,i7 Q \s o = fo, n (^0) l^o — 
fi admits a regular solution on D if and only if 

/ [W y h - f n(tfg] poe- iU °te°dn = 0, Vy e (z s , z N ), 

J PS 

(Hi) the Cauchy boundary value problem A 7 £7q + e~ AUa (df2' , dJ? ) 7 = 0,Uq\z o = fo, 
n({7g ) \s = fi admits a regular solution on D if and only if 
rPN 

/ bPy fi ~ fo n(& v ) - Ti (n)] p dfj, = 0, VyE(z s ,z N ), 

J PS 

and (iv) the Cauchy boundary value problem A 7 .!? ' — 8 (d/2 , dU^)^ — 4 (dJ? ', df/o) 7 = 
0, J?o \s = /o, n (fi'o) \s a — fi admits a regular solution on Do if and only if 

(#„ fi - fo n(!iy) e- 4U ^o _ j 2 (n)l Po dp = 0, Vy e (z s , z N ), 



rPf 
J ps 



where ip y is given in (6), <P y = -^====, T : = ^<2+*dQ_ andT 2 = j;Q-*d (W y + U' Q ) 
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